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ABSTRACT 

Though the bulk of the observed optical flux from the discs of intermediate-redshift lensed 
quasars is formed well outside the region of strong relativistic boosting and light-bending, 
relativistic effects have important influence on microlensing curves. The reason is in the di- 
vergent nature of amplification factors near fold caustics increasingly sensitive to small spatial 
size details. Higher-order disc images produced by strong light bending around the black hole 
may affect the amplification curves, making a contribution of up to several percent near max- 
imum amplification. In accordance with theoretical predictions, some of the observed high- 
amplification events possess fine structure. Here we consider three putative caustic crossing 
events, one by SBS J1520H-530 and two events for individual images of the Einstein's cross 
(QSO J2237H-0305). Using relativistic disc models allows to improve the fits, but the required 
inclinations are high, i > 70°. Such high inclinations apparently contradict the absence of 
any strong absorption that is likely to arise if a disc is observed edge-on through a dust torus. 
Still, the high inclinations are required only for the central parts of the disc, that allows the 
disc itself to be initially tilted by 60. .90° with respect to the black hole and aligned toward the 
black hole equatorial plane near the last stable orbit radius. For SBS J1520H-530, an alterna- 
tive explanation for the observed amplification curve is a superposition of two subsequent fold 
caustic crossings. While relativistic disc models favour black hole masses ^ lO^'^Mo (several 
times higher than the virial estimates) or small Eddington ratios, this model is consistent with 
the observed distribution of galaxies over peculiar velocities only if the black hole mass is 
< 3 X 1O*^M0. 
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1 INTRODUCTION 

The history of gravitational lensing may be traced down 
to the beginning of the X Xth century and beyond (see 
ISchmidt & Wambsgansa ( 120100 and references therein for a histor- 
ical review). First double quas a r imag es produced by strong lenses 
were reported in IWalsh et al.l (Il979h . The first quadruply-lensed 
quasar, QSO J2237+0305, was discovered by|Huchra et al. ( 1985). 
This object was also the first were eff ects of microlens ing by the 
stars of the lensing galaxy were found dlrwin et al.lll989l) . 

The principal difference between strong lensing and mi- 
crolensing is in the angular distance scale set by the Einstein- 
Chwolson radius: 
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Here, A/ is the mass of the lensing object, Dl,s,ls are an- 
gular size distances toward the lens ("L"), the source ("S") and be- 
tween the source and the lens ("LS"), Dls = Ds — Dl x (1 4- 
2l)/(1 + 2s) for a flat Universe (see for example iHogd ( 119991) ). 
Dependence on the lens mass leads to drastically different angular 
scales associated with strong lensing by galaxies and galaxy groups 
(arcseconds) and microlensing by individual stars (microarcsec- 
onds and less). The individual images formed in the latter case 
cannot be resolved by contemporary instrumentation, but their am- 
plification variations may be a valuable tool to resolve the spatial 
structure of the source (Chang & Refsdal 1984) . Accretion discs 
around supermassive black holes at cosmological distances should 
have comparable or somewhat smaller angular sizes at A ~ 2000A. 
Below we will use spatial sizes projected onto the picture plane at 
the distance of the source. They differ from the angular sizes by the 
dimensional factor of Ds '■ 
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The main difference between microlensing by stellar mass ob- 
jects of our Galaxy and microlensing effects accompanying strong 
lensing of quasars is in the optical depth. Angular distances be- 
tween individual lensing bodies in a galaxy scale with the distance 
D as oc D~^, while Einstein-Chwolson radii decrease only as 
oc D""'^. If we define the microlensing optical depth as the to- 
tal solid angle of the Einstein circles of th e stellar population of a 
unit s olid angle of the lensing galaxy (cf. iNitvananda & Ostriken 
( Il984) ). it may be written as: 
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Here, E is the total surface density of the clumped matter 
in the lensing galaxy. Smoothly distributed component (primar- 
ily, unclumped dark matter) contributes only to strong lensing, un- 
less it creates a strong shear that destroys point-lens degeneracy 
JNitvananda & Ostriken 1 19841) . Optical depth becomes consider- 
ably large (r > 0.1) forlensed quasar^at gigaparsec distances. 

Difference in optical depth makes microlensing effects in dis- 
tant lensing galaxies qualitatively different from the rare single- and 
double-lens events in our Galaxy. As a background source moves 
with respect to a single point-l ike lens, it has z ero probability to 
undergo infinite amplification jPaczvnskil 1 1 986 b ). In this regime, 
amplification becomes sensitive to the size and the structure of the 
object only if its angular distance from the lens centre becomes 
comparably small. On the other hand, strong (r > 0.1. .0.5, de- 
pending on the underlying shear) microlensing by a population of 
point masses creates a network of fold caustics jPaczvnskill 1 986al) 
where amplification is divergent and beh aves as oc d~^'^, where d 
is the angular distance toward the fold JChang & Refsdallll984l) . 
In the case of quasar microlensing, every particular image tra- 
verses some fold caustic at a probability of about unity on a several 
years' time scale (for the case of QSO J2237-I-0305, the relative 
fre quency of high-am plification events was estimated as ~ lyr~^ 
by Iwitt et alj ( 1 1993b . other objects have smaller relative proper 
motions). A comprehensive review on microlensing was made by 
[Wambsganss (2006). 

At present, there are at least three different approaches to 
quasar microlensing that may be used to probe the structure of 
supermassive black hole accretion discs. First is to gather statis- 
tics on image amplification and to compare single-epoch anoma- 
lous flu xes with the predic tions of accretion di sc models, as it wa s 
donebvlPoolev et aljj200 7): Bate et al. (2008): lFlovd et alj ( l2009l) : 



iBlackbume et al.l ( 1201 ll) : I Jimenez- Vicente et alj ( l2012h . In the first 
paper, a huge, about a factor of 10.. 100, inconsistency was found 
between the sizes of accretion discs estimated by microlensing 
methods and predictions of the standard accretion disc theory. Par- 
tially, this inconsistency may be attributed to the crude mass esti- 
mates applying broad-band magnitudes with some assumptions on 
bolometric correctio n and accretion efficiency. On the other hand, 
[Morgan et al.l ( 120101) use more accurate mass estimates based on 
the widths of broad emission lines and yield much better consis- 
tency. Indeed, among the four objects common for the samples 
used in these two studies, only for one (QSO J2237-I-0305) the 
masses determined by the two methods are consistent within the 



^ Here we do not distinguish between radio-loud quasars and radio-quiet 
"quasi-stellar objects" (QSO) and refer to both object types as "quasars" or 
QSO. 



uncertainties. For two objects (PG 1115-1-080 and RXJ 1131-1231), 
photometry-based mass estimates are about four times lower, while 
for the least massive one, SDSS 0924-1-0219, masses differ by an 
order of magn i tude. However, the accretion disc sizes measured by 
[Morgan et al.l ( 120101) are still several times larger than expected. 
Authors suppose that virial mass estimates may still be systemati- 
cally lower by a factor of ~ 3. 

The analysis method u sed bv lMorgan et alj ( l20ld) was intro- 
duced by iKochaneM (20040 and may be characterised as extensive 
light curve fitting. A number of artificial amplification maps is gen- 
erated, and the observational light curve is compared with numer- 
ous simulated light curves with random parameter s. This technique 
was ap plie d in a large number of studies such as lEigenbrod et al.l 
( 120081) and lHainline et alj ( l2012l) . Since the technique requires mul- 
tiple observational points, it is perfect for Einstein's cross and other 
targets of monitoring programmes. 

The Monte Carlo microlensing analysis is rather resource- 
consuming. As an alternative to computationally-extensive meth- 
ods, it is reasonable to study individual high-amplification events 
that are most sensitive to the spatial properties of the source. 
Primarily, high amplification events are associated to fold caus- 
tic crossings, when a pair of new microlensing images ap- 
pears or disappears, and the point-source amplification diverges. 
Caustic crossings by standard discs are good mod els for some 
of the observed am plification events (Gil-Merin o et al.l 120061 ; 
iKoptelova et alj|2007l) . It may be shown that point- source ampli- 
fication during a caustic crossing is divergent for a disc with a 
small inner radius rin as /i oc r -^ , hence brightness maxima in 
microlensing curves are best fit f or resolving the innermost parts 
of the disc 1 Agol & KroliMll999. ). Indeed, as it was shown by 
Ijaroszvnski et al.l ( 119921) , lensing curves differ for accretion discs 
of different inclinations and Kerr parameters of the accretor. Am- 
plification curves primarily differ near their maxima. 

Caustic crossing events allow to study the structure of the in- 
nermost parts of accretion discs where general relativity effects are 
important. Two processes are expected to influence considerably 
the observed light curves: light bending and Doppler boosting (due 
to matter motion in the disc as well as due to frame dragging). Both 
make brightness distributions asymmetric and strongly dependent 
on the inclination angle. In this paper we aim on estimating the in- 
fluence of relativistic effects on the amplification curves created by 
straight caustic crossing events. We also apply the results of our 
calculations to three high amplification events and show that some 
of their features are probably connected to relativistic effects. 

The paper is organized as follows: first we describe the 
archival observational data we use. Simulation technique applying 
Kerr geodesic calculation software is considered in section [3] Fit- 
ting results for the three putative caustic crossings are given in sec- 
tion |4]and discussed in section |5] 



2 OBSERVATIONAL DATA 

The current number of lensed qu asars where microlensing effects 
were found is about several tens (Pooley et aLluOOTI : [Morgan et alj 
[2010; Jimenez- Vicente et al. 2012). However, long homogeneous 
observational series exist for few objects. The best studied among 
these is the Einstein cross (QSO J2237-I-0305) that is a subject 
of extensive photome tric monit oring programs such as OGLE- II 
( IWozniaketal] |2000^) and -III JUdalski et al] |200d^. We analyse 
two high-amplification events that took place for two different im- 
ages in the years 1999 and 2000. We also found a considerable 
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amount of archival data on another object, SBS J 1520+530, and 
interpret variations of the relative image amplifications as a man- 
ifestation of microlensing amplification. A brief summary of the 
observational data we use in this work is given in table[T] 

2.1 SBS J1520+530 

The object is doubly ima ged, with the avera ge flux ratio for the two 
images A and B about 2 i Burud et al.l2002!) . The lens is a relatively 
distant (z ~ 0.7) late-type elliptical with the estimated velocity dis- 
persion of (J ~ 200kni s ~ ^ , steep mass profile, central convergence 
of K :: ::: 0.5 and possible signatures of interaction with its environ- 
ment JAuger et al.ll2008h . B is about three times closer to the lens 
centre, hence we expect moderate microlensing optical depth for B 
and small (r < 0.1) for A. Hence we are inclined to interpret the 
observed flux ratio variations as microlensing of the B image. 

We used the three following sources of reduced photometric 
data (R-band magnitudes for both images) covering a time span of 
about seven years (see figure[T): 

• 58 data points obtained with the Nordic Optical Telescope 
(NOT) and published by Burud et al. (2002) 

• 60 observations with the 1.5-m Russian-Turkish Telescope 
RTT jKhamitov et alj|200a) . kindly provided by I. Bikmaev and 
I. Khamitov 

• 123 observations with the 1.5m AZT-22 telescope in Maid- 
anak. These data wer e described, analysed and published by 
iGavnuUina et alj ( |2005|) . 

All the photometry was performed in the standard optical R 
band that for the source redshift of 1.855 ( iBarkhouse & Halll200lh 
corresponds to A ~ 2000..2500A in the source reference frame. 

Relatively large amount of data allows to improve the esti- 
mate of delay time. Cross-correlating the time series for the two 
images we find a broad peak at At = 127.6 ± 2.0 days that we in- 
terpret as the delay between the two images. Within the uncertain- 
ties (that are hereafter calculated for 90% significanc e level), this 
value is consistent with the 130±3 day delay found bv lBurud et al.l 
( 12002) and Gaynullin a"et alj (|2005l) and with the 128-day estimate 
bv lKhamitov et alj (2006). This value of At is subsequently used 
to shift the series of A and B fluxes and estimate flux ratios. 

Merged time series consists of 241 unevenly-sampled R-band 
observational points for each image (see figure [T). To minimize 
information losses but to avoid unjustified interpolation (on pe- 
riods of time longer than tens of days), we construct flux ratios 
by linearly interpolating the magnitudes between the observational 
points. The resulting flux ratio series consists of points of two 
kinds: (i) B image fluxes shifted backwards by At and divided 
by the interpolated values of A image flux and (ii) interpolated B 
image fluxes divided by the A image fluxes shifted forward by At. 
We did not interpolate over time gaps longer than 45d (about 16d in 
the frame co-moving with the object), therefore the resulting time 
series contains only 360 points, not entirely independent (see also 
section l4.U . 

Flux ratio curve (see figure|2) demonstrates a maximum at the 
Julian date JD ~ 2451500 and another one at JD ~ 2451800, 
that form together an about two-year long structure similar to stan- 
dard disc amplification curves (see next section). For the remain- 
ing two-thirds of the curve, B/A ratio is stable and close to 0.47 
with standard deviation around four percent. From the resulting 
time series, we excluded the final portion of the curve between 
JD = 2453000.. 2453600 where the data are sparse and the ob- 
served variability is difficult to relate to the high-amplification 



event in the earlier part of the data (see figure 1101 . The resulting 
amplification curve consists of 253 data points. 

2.2 QSO J2237+0305 

We use R- and V-band photometric data on the images A and C of 
QSO J2237+0305. The data were taken from two sources: 



Huch ra's lens monitoring program 
|200C|) , the data are available at 



• OGLE-II 
dWozniak et al.l 

http: //ogle. astrouw. edu .pl/cont/4_niain/len/huchra/huchra_ 
and contain V-band ma gnitudes for the four im ages 

• GLITP archive (Alcald e et all |2002|) : V- and R-band 
magnitudes for all the images, reduced data are available at 

|http ://wela.astro.ulg.ac .be/themes/extragal/gravlens/bibd 

The light curves are shown in figure [3] The data sample is 
similar to that used bv lKoptelova et alj ( 120071) . We use GLITP data 
reduced by two different techniques: ISIS photometry (the magni- 
tudes are available at the web page given above) and PSF-fitting 
(the magnitudes are meant to be available for download but the 
link is broken, hence we asked Elena Shimanovskaya who has 
kindly provided us with these magnitudes). T he two reduction tech- 
niques are compared in lAlcalde et al.l (120021) . Magnitudes reduced 
by different techniques are generally consistent but deviate consid- 
erably (by about 0.™03) near the maximum of the image A high- 
amplification event. Unfortunately, the maximum interferes with 
a period of bad visibility of the object (January/February 2000). 
OGLE points show unreasonably large scatter near JD=2451500. 
Besides, consistency between OGLE and GLITP data during this 
period is poor and insufficient for our purposes (see figure |3] left 
panel). Better but still significant deviations between the GLITP 
magnitudes reduced with different methods (we show only ISIS- 
reduced data to prevent the figure from overcrowding; the two 
GLITP light curves are shown together in figure 1121 . Therefore 
we exclude OGLE data from subsequent analysis for the image A 
event. 

GLITP data during the maximum of image A brightness 
demonstrate a double-peaked structure that we are inclined to in- 
terpret as a signature of the inner disc structure. Overall behaviour 
near Julian dates 2451520. .2451540 suggests that the dip is real 
and has an amplitude of about 0.™02. 

For the C image event, most of GLITP observational points 
are far away from the high amplification event maximum, and we 
finally fit only the V-band OGLE data. 

There is strong evidence that the individual time delays for the 
ima ges of QSO J2237+030 5 are very small, such as several days 
(see Koptelova et alj (200g) and references therein). We can make 
use of the relative proximity of the lens that makes the microlensing 
variability of individual images much faster and stronger than the 
intrinsic variability of the source. The studied high amplification 
events last less than a year, and variability of individual images does 
not show any detectable correlation on these time scales. Einstein 
cross is known for its lack of variability, the maximal gradients 
of intrinsic variability are considerably smalle r than microlensing 
trends (see for example lEigenbrod et alj(l2008h ). 



3 AMPLIFICATION CURVE SIMULATION TECHNIQUE 

3.1 Basic assumptions 

All the three high-amplification events are good candidates for 
caustic crossings. We model them by convolving a straight caustic 
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Table 1. Basic information about tlie objects and the observational amplification curves used for analysis. For the QSO J2237+0305 image A event, the 
numbers of GLITP data points for two reduction techniques, ISIS and PSF, are given separated by a slash. 





SBS J1520+530 


QSO J2237+0305(A) QSO J2237+0305(C) 


source redshift 
lens redshift 


1.855 
0.72 


1.695 
0.039 


time span (V), JD-2450000 
number of points (V) 


- 


1400.. 1650 1200.. 1650 
53/52 83 



time span (R), JD-2450000 1200.. 3000 1450.. 1510 

number of points (R) 253 51/49 



2000 



2002 



year 



2004 



2006 



8.10 

8.20 

8.30 

8.40 

8.50 

8.80 
8.90 
9.00 
9.10 
9.20 
9.30 
9.40 
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Figure 1. R-band light curves for A (upper panel) and B (lower panel) images of SBS J1520+530. Filled circles, asterisks and diamonds represent NOT, RTT 
and Maidanak data (see text, section[2]TJ. 



crossing curve with a standard disc profile integrated over one di- 
mension. Simplified flat and non-relativistic (section [T2t and gen- 
eral relativistic (section [53) disc models were used. In both cases 
we consider a multi-blackbody disc with monochromatic intensity 
obeying the Planck law with the temperature determined by the 
standard accretion disc model: 



/ex 



^Ri^ 



expr(r/rd)3/4(l-yr,„/r)^ 



(4) 



where r is the radial coordinate, and rd and ri„ are, coiTe- 
spondingly, the radial scale of the disc and its iimermost stable orbit 
radius defined as: 



rd = 



A 



he 1 + z 
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, ,GM 
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(5) 



(6) 



Here, k, a and h are Boltzmann, Stephan-Boltzmann and 
Planck constants, c is the speed of light, / and rj are Edding- 
ton ratio and overall accretion efficiency, x is Thomson opacity. 
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Figure 2. Ratios of delay-corrected fluxes of images A and B of the lensed QSO SBS J1520+530. Tliin liorizontal line is the mean ratio value, thick horizontal 
line with ticks marks the time span used for our analysis. 
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Figure 3. Raw light curves of the images A (left) and C (right) of QSO J2237-I-0305 used in the cunent work. GLITP(ISIS) data ai'e shown by filled circles, 
OGLE data points by asterisks. Date intereals used for fitting are restricted by vertical lines. 



K ~ 0.35cin^ g^^ for Solar composition. Dimensionless inner- 
most stable circular orbit radius 2;isco(o), as well as efficiency 
r] = r]{a), may be found analytically a s functions of dimensionless 
Kerr parameter a dSardeen et al .11 19721) . 

Expression (|4} was derived for monochromatic intensity (no 
matter I,j or Ix) but is to a high accuracy valid even for broad-band 
photometry if most of the radiation comes in the form of a smooth 
continuum. Mean effective wavelength varies by less than 2% for 



the spectral slope varying from p = — ltop= 1, where F^ cc v^ . 
The ab ove definition of r^ i s identical to the characteristic radius 
used bv lMorgan et alJOOlOf) . Note however that it is generally suf- 
ficiently smaller (by a factor of ~ 2.44 for a standard disc with no 
inner edge) than the half-light radius. We will see below that for the 
shape of amplification curve, the quantity of principal importance 
is the ratio of the two radial scales given by (|5) and (|6): 
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Figure 4. A qualitative scheme showing the two coordinate systems in the 
picture plane. Shaded is the area behind the caustic where amplification is 
higher. Disc is inclined with respect to the observer, and its closer side is 
shown by a thicker curve. 
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For fixed X, the differences between amplification curves 
are important only for high inclinations and are connected to 
relativistic effects, namely Doppler boosting and light bending. 
Variations of the Kerr parameter may change X by a factor of 
msQ<i(xisco)/ Tain(xisco) = 9. Inversely, one may reasonably 
estimate X by fitting the amplification curve and arrive to a tightly 
correlated pair of uncertain a and M. Having reasonable mass es- 
timates, one may thus make an estimate for a, and vice versa. 

Below, we fix the Eddington ratio to I = 0.25. The black hole 
masses are estimated as M = 8.8 x 10* M© for S BS J1520+530 
and M = 9 x lO^M© for QSO J2237+0305 bv [Morgan et"aLl 
( 120100 . These mass estimates were obtained using emission line 
profiles (see references given by Mor gan et al.lcOlOl) ) and virial re- 
lations (jVestergaard & Peterson 200q). Accuracy of these estimates 
is low, about 0.3dex. Our analysis yields generally higher masses 
poorly consistent with the virial estimates (see below section|4j. In 
fitting the observational data, we first find X, inclination i and po- 
sitional angle ^ and then optimize for the mass and Kerr parameter, 
fixing X. 

Caustic is considered a straight line defined by condition 
y ~ '"eff X (t — to) at the given time t. Let x and y be the co- 
ordinates in the picture frame along and across the caustic, respec- 
tively. The origin of this coordinate system coincides with the black 
hole centre. In the general case, the disc is inclined, and it is con- 
venient to use the coordinate system (a, /3) defined by the major 
and t he minor axes o f the d isc projection upon the picture frame 
(as in lOexter & Agoll ( 120090 ). The two coordinate systems in the 
picture frame are connected by rotation by an angle of ^p that has 
the meaning of relative positional angle of the normal to the disc 
with respect to the normal to the caustic (see a sketch in figure |4ll. 
Without any loss of generality we assume that the relative motion 
of the source and lens is normal to the caustic itself. Integration 
over one direction is convenient to consider as projection upon the 
orthogonal direction. 

Effective velocity v^ff is calculated as the relative proper mo- 
tion multiplied by the angular size distance Ds toward the object. 
Our analysis is sensitive only to its component perpendicular to the 



caustic direction. The velocity is connected to the peculiar spatial 
velocities of the QSO (vs) and the lensing galaxy (vl) in the fol- 
lowing way: 



"e// = 



VS n 



Vl 



n Dc 



+ 



Vo • n Dls 



1 + ZS 1 + ZlDl l + ZL D 



(8) 



Third term describes the peculiar motion of the observer, 
Vo ~ 370km s^^ is known from the di pole constituent of cosm o- 
logical CMB variations (see for example iLineweaver et al.l ( 119960 ). 
Unit vector n is normal to the caustic and lies in the picture 
plane. Reasonable limit for peculiar motion velocities of indi- 
vidual galaxies is « < Vmax — 2000kins~^ that corresponds 
to abou t 2.. 3 rms values for the radi al peculiar velocity distri- 
bution JRavchaudhurv & Saslawlll99q ). This limit may be con- 
verted to the condition for the possible effective transverse veloc- 
ity Veff < V„,aa:^/{1 + Zs^ + (1 + Zl)-'' X (Ds/Dl)^. For 

SBS J 1520+530, vs and vl contributions are of the same order 
because Ds/Dl — 1.2. In this case, v^ff < 1600km s~^. For 
QSO J2237+0305, the lens is about ten times closer (Ds/Dl ~ 
ll)andi;e// < 20 000kms"\ 



(7) 



)(«) 



3.2 Simplified standard disc 

Besides the more sophisticated model described below, we use a 
simple disc approximation ignoring all relativistic effects. In this 
approximation and for the straight caustic case, inclination does 
not affect the observed shape of the amplification curve. This is 
valid for a thin disc of arbitrary inclination as long as the disc is 
thin enough (h/R ^ cos i, where i is inclination). 

We are interested in intensities integrated over one dimension. 
Here, we provide the general form for the integral over one direc- 
tion (we make a substitute t — x/y): 

h{y) = J I{x,y)dxcx 



vi: 



3/S 



where: 



K = K{i,'ip) =cos^Tp + 



sin ip 



■ 2 , , cos^ i) sin^ V cos^ tp 4 . 
J — J[i,ip) = sm %p H — 777^ — 77 — tan i, 



K[i,^) 



1 r^ 



1 



J y2 1 + i2 



1/4 







T^in — 



iiy^r^n/^fj 



(9) 

(10) 
(11) 
(12) 

(13) 



If the influence of the inner radius is negligible, the one- 
dimensional intensity profile is identical to the face-on disc inten- 
sity profile with the radial scale of r^ = rd/vJ- Even if the cor- 
rection factor / is taken into account, profile shapes do not depend 
on the angles ijj and i, because the integral (|9) contains the angles 
and coordinates only in combinations yyJ jrd and yvJ /tin that 
affects only the stretch factor of the curve for given X. 

The influence of the inner disc edge is included by the single 
parameter X = rd/vin (see previous subsection) determining the 
shape of the amplification curve. In figure |5]the dependence of the 
amplification curve on this parameter is shown for a fixed black 
hole mass. Kerr parameter is varied from -0.99 to 0.99. 
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3.3 Kerr black hole disc of arbitrary inclination 

Photons move along null geodesies that we calculate using 
geokerr code JDexter & Agolll2009r) . Disc was ascribed a con- 
stant relative thickness of h/r — 0.01 that does not affect the re- 
sults much as long as cos i > h/r. The maximal inclination value 
we use is i = 88° that implies cos i ~ 0.03. We linearly interpolate 
the radial coordinate to the point where cos i — ±h/r to determine 
the intrinsic intensity of the locus in the disc corresponding to cur- 
rent values of a and /?. For this point, using the conserved angular 
momentum k^ and the energy at infinity kt of the photon, we con- 
nect the observed photon energy (equal to kt) to its energy in the 
frame co-rotating with the disc. Doppler- factor is calculated as: 



Vpba 



kt 
kiU' 



(1 + fll) 



(14) 



Here, u* = l/\/—{gtt + 2gtipfl + g^^Jl^) is the null com- 
ponent of the disc rotation four- velocity (where metric signature is 
-l~l~l-), I = k^/kt = asini is the net angular momentum of the 
photon. For co-rotating discs. 



n 



r3/2 + , 



(15) 



in c^ /GM units, while the radius r is in GM/c . Counter- 
rotating discs may be also considered in this formalism by changing 
the sign of a. Hereafter we will use negative a values to describe 
the case of counter-rotation. 



Observed intensity is integrated over a fixed (observer-frame) 
wavelength range and is therefore approximately proportional to 
the monochromatic intensity I at the central observer-frame fre- 
quency Uoba- For a local blackbody-like spectrum (we consider /„ 
for convenience): 



oc exp i [r/r-d) 1 - \fT~Jr- 



1/4 



(16) 



Doppler boost thus contributes only through direct frequency 
shift. At large distances (r 2> Tin and r ^ r^) its effect on the 
local intensity is still important, primarily due to the rapid fall- 
off of the Planck law with frequency. For inclined discs, intensity 
distribution along the major axis is asymmetric even at large dis- 
tances from the black hole. Up to the two leading teims in 1/r, 
1/(5 ~ 1 + sin(i)/y'r, and relative intensity difference between 
the approaching (/"*") and receding (/~) sides of an inclined disc 
is: 

7+ -7 



7+ + 7- 



(r/rd 



3/4 sm(zj 



1/4 



(17) 



Above estimate was made in the assumption that the asymme- 
try is small. The actual value of asymmetry is evidently limited by 
the value of max ((7"*" — 7~)/(7+ + 7~)) = 1 when one side of 
the disc is much brighter than the other. The approaching side of an 
accretion disc is generally about two times brighter than the reced- 
ing at r ~ 10..100GA7/c^. Weighed disc centre at large distances 
is effectively shifted by Ar ~ sini x r^ ". 
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Figure 5. Amplification curves for a simplified disc with a fixed mass 
(M = lO^M© at comoving A = 2300A) of the black hole and variable 
Kerr parameter. 
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3.4 Calculation of amplification curves 

For calculating the shapes of null geodesies, we use the public code 
oflDexter & Agol (2009). Then we take into account Doppler shifts 
and Doppler boosts in the way described above. For given Kerr pa- 
rameter a, inclination i and relative positional angle ijj, intensity is 
calculated on a rectangular grid log-uniform in x and y coordinates. 
Then the two-dimensional brightness distribution is integrated over 
one dimension: 



h(y) = 



I{x,y)dx 



Then, the microlensing amplification curve for the straight 
causti c case may be calculat ed in the following way (see for ex- 
ample |Ko2telova_et_sLl (20071) who used a similar approach): 

m ^ ^ J hjy) X Ay-'^' X e{Ay)dy , 
^(f ) = ^1 V Co X ^^ ,. ^ , ^ , h^to (18) 



Ihiy)dy 

Here, Ay — y — v{t — to), and 0{x) is Heaviside function: 
Q{x > 0) = 1 and8(x ^ 0) = 0. Spatial scale factor (o is chosen 
according to lWitt et al] ( Il993l) : 



Co = 



AGMr. 



il-K, 



.DlsDs 

' Dr 



= rEi 



,{Mq) X ^/^^^^ (19) 



Here, Kc is the continuous contribution to the total con- 
vergence. If we adopt Kc = 0.5, for SBS J1520-I-530 and 
QSO J2237-I-0305, the spatial scale is Co i; 3 x lO^'' and 1.3 x 
lO^^cm, respectively. 

Effective transverse velocity v and the fold-crossing epoch to 
are considered as free parameters. For every model (every given 
one-dimensional brightness distribution), v and io are optimized 
for by minimising x^ for the observational amplification curve. For 
the A image of QSO J2237-I-0305, two amplification curves were 
used for the two filters, R and V, and we fit both simultaneously. 
Our fitting has two additional parameters, /io and ni, that we cal- 
culate using linear regression for each iteration of the optimisation 
process. In the case of QSO J2237-I-0305, two /io,i pairs were used 
for the two photometric bands. The ratio of the two coefficients, 
Hi/fiQ, has the meaning of cau stic strength an d depends only on 
the mass distribution in the lens dKavser & Wit t 1989). 

Evidently, for the shape of the caustic-crossing amplification 
curve, the effects of general relativity are of primary importance. In 
figiire|6] we show four light curves emerging from caustic crossings 
by a disc inclined by 60° around aa — 0.6 black hole, traversing a 
straight fold caustic in four directions: along the minor axis (tp — 
and 180°) and along the major (tp — 90 and 270°). Effective trans- 
verse velocity and black hole mass used for the figures |6] and |7] 
were 5000km s~^ and 10^ M©, respectively. Due to their asymme- 
try, inclined discs produce amplification curves strongly dependent 
on the relative positional angle. Figure |6]demonstrates all the main 
complications arising from relativistic effects. If the brighter part of 
the disc is amplified while its dimmer part is still unaffected by the 
caustic, amplification curve may become two-peaked with the rel- 
ative intensity of the two peaks strongly dependent on inclination, 
positional angle and black hole rotation parameter (figure|7). An ac- 
cretion disc surrounding a rapidly rotating Kerr black hole observed 
at high inclination has a bright compact "hot spot" that influences 
strongly the amplification curve. Relativistic effects are capable, in 
particular, to qualitatively explain the observed fine structure of the 
amplification events under consideration. 

Test runs reveal a small but considerable contribution from the 
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Figure 6. One-dimensional brightness distributions (upper panel) and am- 
plification curves (lower panel) for a disc inclined by 60° around a a = 0.6 
black hole. Dashed curve is the simplified disc model, black and grey curves 
correspond to caustic transition along the minor and major axes of the in- 
clined disc (solid and dotted curves differ in propagation direction). 



photons belonging to the higher-or der images produced by e xtreme 
light bending near the black hole JBeckwith & Donell2005n . Their 
contribution to the total flux may be as high as several percent, 
depending on the in clination and on th e Kerr parameter. The fig- 
ure of 10% given bv lBeckwith & Done! ( 1200 5) is an over-estimate 
for an optically thick disc because at high inclinations and high 
Kerr parameters considerable part of the black hole is covered by 
the disc. Simulated black hole shadows with one higher-order im- 
age visible are shown in figure [8] Higher orders have significantly 
smaller fluxes due to gradually decreasing solid angle. Flux of n-th 
image loop decreases proportionally to its radial size Ar (inten- 
sity is approximately conserved). During caustic crossing, the ad- 
ditional amplification factor is proportional to Ar~^'^ , where the 
observed projected width of the approximately annular image is 
oc Ar. The resulting contribution of the n-th order image scales as 
oc Ar^ " oc exp(— TTTi/b), where h ~ 2.7 (in the so-called strong 
deflection limit, see iBozzal J20ICI) ). Higher-order images are thus 
exponentially damped. 

An example of higher-order contribution is given in figure |9] 
For high inclinations, it is generally about a couple percent, but 
are less than one percent for i ~ and a > 0. For counter-rotating 
black holes, the size of the inner hole is the largest and the contribu- 
tion of the secondary images reaches 2.. 3 percent for characteristic 
QSO black hole masses ~ 10^ Mq. The effect is also important for 
larger black holes and at shorter wavelengths. 
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Figure 7. Amplification cui'ves for an accretion disc around a black hole with a = 0.6 inclined by i = 60° (upper panel) and i = 85° (lower panel) for 
variable positional angle ^. Every horizontal sUce corresponds to one curve. Scale is linear, maximal shade con'esponds to maximal brightness. 



4 RESULTS 

Fitting results with the simplified face-on disc and fully relativistic 
disc models are given in tables |2] and [5] respectively. 



4.1 SBS J1520+530 

Fitting with the plain face-on disc model allows to approximately 
recover the size of the inner hole. Number of degrees of freedom 
is 253 — 5 = 248 here ("optimistic" number, see below). We fixed 
the mass of the black hole to 3 x 10^ M0 and localized an abso- 
lute minimum near a = —0.4. Corresponding X value is 1..2. 
The minimum is consistent with physical — 1 < a < 1 only 
for the black hole masses in the range (1.6. .3.5) x lO^M©. The 
best-fit amplification curve is also shown in figuresfTOlandllllbv a 
dashed line. Here we limited the effective velocity by the value of 
Veff,max = 2000km s~^ . If We rclax thls limitation, the fit may be 



improved (up to x^/DOF ~ 190) by increasing the mass of the 
black hole and the transverse velocity. 

A total number of 708 fully relativistic models was calculated 
covering the possible ranges of inclination < i < 90°, relative 
position angle < ■(/) < 360° and Kerr parameter for co-rotating 
(1 > a > 0) and counter-rotating (— 1 < a < 0) cases. For every 
model, we make an optimization run finding the best-fit v^ff and to 
and two amplification parameters /io and fii . Thus the model one- 
dimensional brightness distribution fixes only the shape of the am- 
plification curve that may be then stretched and shifted in time (to 
and Veff) and in amplification factors (/io,i). For Kerr disc models, 
we first fix the mass of the black hole to 0.88 x lO^M© and then. 



after finding the global minimum at a = —0.9 



+0.4 



I ±5° 



and Tp — 340 ± 4°, fix X, i and 'ip and vary the mass and rotation 
parameter. 

In tables |2] and [3] we give the best fit parameters for the am- 
plification curve of SBS J1520+530 fitted with the simplified disc 
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a=0 



a=0.9 





Figure 8. Black hole shadows for i = 80° and two Kerr parameters, a = and 0.9. Secondary image is clearly seen for Schwarzschild black hole and is 
visible as an inclined nearly straight line at the approaching (left) side in the other case. A bar below the black hole has the length equal to the innermost stable 
orbit radius (6 and ~2.3 in GMji? units, respectively). 



model and with a Kerr disc with the best-fit parameters. Fully 
relativistic disc provides much better fit (see below this section). 
However, the observed details in the amplification curve are even 
sharper and more profound. This is a possible signature of a non- 
trivial structure of the inner parts of the disc resulting from either 
disc tilts and warps (see below sect ion [54t or from additional en- 



ergy input from black hole rotation dAgol & KrolikluOOOf) . 

The apparently low values of x^ should be addressed sepa- 
rately. The procedure that we used to reproduce the amplification 
curve (see section |2ll minimizes the loss of observational informa- 
tion but generates two sets of flux ratios (A image fluxes interpo- 
lated over the observational epochs for B and vice versci) that are 
not entirely independent. Therefore, the effective number of de- 



grees of freedom (DOF) is not the number of flux ratio data points 
CNf =253) minus the number of parameters {Np — 5 for sim- 
plified disc and Np = 8 for Kerr metric) but less. Therefore in the 
tables we give two estimates for the number of DOF: one optimistic 
DOFi = iVF-A''p, and one pessimistic DOF2 = Nf/2-Np ~ 
120. The real value is somewhere in between. For the simplified 
disc model, the number of parameters is less by 3 (i and ^ lacking, 
and X instead of a and M). 

For both DOF estimates, the difference in x^ is significant: F- 
test gives the probability of 7 x 10^^° (optimistic) and 4 x 10~^ 
(pessimistic DOF estimate) for the more sophisticated model to 
provide a better fit accidentally. The measured best-fit parameters, 
primarily the high inclination, are suspicious but not entirely un- 



Figure 9. Relative difference between the microlensing curves calculated 
in the single-image assumption and taking into account all the resolved sec- 
ondary loops. We set M = IQIOMq, a = 0, J = 80° and V = 
(solid line) and 90° (dotted line). Effective transverse velocity is fe// = 
1000km s~^, vertical line marks the instance of caustic crossing by the 
centre of the black hole. 
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Figure 10. R-band B/A flux ratio for SBS J1520+530 fitted with a Scliwarzscliild face-on "plain disc" witliout relativistic effects (dashed curve) and a Kerr 
multi-blackbody disc of arbitrary inclination and spin (solid). The parameters of the best-fit Kerr disc model are: a = —0.5, M = 2 X IO^Mq, i = 80°, 
^p = 330° . Lower panel shows residuals with respect to the Kerr disc model. 



physical. We propose that the real disc should be inclined with 
respect to the black hole plane and its inner parts should be thus 
strongly warped and distorted. In more detail we discuss this issue 
in section[5!4l 



An alternative to a single fold caustic amplifying a relativis- 
tic strongly-inclined disc is a more complex structure of several 
caustics or cusps. We find excellent agreement fitting the ampli- 
fication curve with a simplified disc traversing two fold caustics 
sequentially in one direction at one velocity (figure 111) . Best-fit 
parameters are, for fixed a = 0.2 and M = 8.8 x 10*Mq: 
Veff = (1.35 ± 0.22) X 10''kms"\ toi = 2451440 ± 12d for 
the first fold and to2 = 2451810 ± 30d for the second, the ratio 
of caustic strengths I = 0.84 ± 0.06, x^ = 160/244(118). The 
strengths of the two caustics are about 0.7 and 0.6, the weaker one 
is delayed by At — 367 ± 8 days. The worst problem for this inter- 
pretation of the observational data is the high required transverse 
velocity. Spatial size of the disc is oc M'^'^, and the mass should 
be sufficiently decreased to match the required transverse velocity 
to the observed peculiar velocity distribution for galaxies. Setting 
M = 3 X 10* M© allows to obtain a reasonable fit (x^ ^ 162) for 
a relatively small transverse velocity v^ff ~ 2 x lO'^kms^^. 



4.2 QSO J2237-I-0305 

Simplified disc fits appear even worse for the case of 
QSO J2237-I-0305. Qualitatively, the situation is understandable: 
plain disc symmetric brightness distribution is unable to reproduce 
the two-peak structure of image A event as well as the narrow peak 
feature of the image C event. Both are relatively well reproduced if 
the inclination is high (i > 70°). 

Running more than 2000 models, we find a reasonable fit for 
the A image amplification curve (figure [T2t . For fitting, we use 
GLITP data reduced by the two standard techniques separately. 
PSF-fitting data have larger (but probably more reliable) statistical 
errors hence the minimal x^ is considerably smaller. Higher rota- 
tion parameters (up to ~ 0.9) and higher masses (up to ~ IO^^Mq) 
produce apparently worse fits but are still allowed if one increases 
the observational uncertainties by a factor of several. Inclination i 
and positional angle ^p are well constrained, because the notable 
two-peak structure requires certain disc orientation (see figure |7]l. 
At the moment we are satisfied with the qualitative agreement be- 
cause the data definitely have some additional error sources pos- 
sibly connected to the bad visibility of the object near the ampli- 
fication curve maximum. Amplification curve shape is similar for 
both filters, even the two maxima seen in the R-band curve are also 
visible in the V-band. 

It is impossible to explain the observed amplification curve 
of the C image without a relatively large a > 0.95. The rapid 
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Table 2. Results of amplification curve fitting with the simplified accretion disc model. Eifective velocity is normalised by the disc radial scale (R-band for 
SBS J1520+530, V-band for QSO J2237+0305) in lOi^cm units. 



r-d/rtn 



^e///»-d,i5.km£ 



to, JD-2450000d 



Ml /mo 



xVdof 



SBS J1520+530 



1.5 



+0.6 



740 



+80 
-60 



1676.3 



+4.7 
-0.1 



1.05±0.17 



209/246(120) 



QSO J2237+0305 (A): 
ISIS 

PSF 



2.0 ± 0.03 (R) 
1.6±0.02(V) 



2.12 



+0.1 



(V) 



-2900±100 
-2920±120 



1460±5 
1455.0±2 



2.4±0.3(R) 
3.9 ± 0.4(V) 



4.5! 
6.0 



+0.09 



°(R) 

(V) 



560/99 
122/96 



QSO J2237+0305 (C) 



>4.4 



6100+:? 



+400 



300 



1389.0l55;^ 



0.66±0.02 



183/77 



Table 3. Results of amplification curve fitting with an inclined Kerr disc model. 



rd/rin 



M 
IO^Mq 



deg 






to Ml /mo 

JD-2450000d 



X^/DOF 



SBS J1520+530: 

2.15±0.05 


-0.6.. 0.2 


1.7. .3.4 


80±10 


335±5 


1500..2400 


1678±5 


0.61 ±0.05 


176/244(118) 


QSO J2237+0305: 

A image, ISIS 
1.7±0.1 (R) 0.2+[;-^ 
1.4±0.1 (V) 


7+2 


71+15 


96±5 


-{3±1) X 10* 


1509.9 ±1.4 


0.91±0.06(R) 
0.92±0.03 (V) 


333/96 


A image, PSF 

i-7to;o2 (R) 

1 -410.02 (V) 


-0.3. .0.8 


5t^ 


™l^S 


□0+30 


-(1.5.. 3.0) X lO"* 


1512^^ 


0.96±0.07 (R) 
1.2±0.2(V) 


98/93 



C image 

2.1+«- 



0.97±0.02 



10^ 



85±5 53±10 (2.0±0.2) X lO'' 



1309±4 



2.09±0.05 



88/75 



rise of the curve in figure [13] can not be reproduced unless the 
accretion disc has a bright spot naturally explained by high rota- 
tion parameter, high mass models. The predicted mass is very high, 
M > 7 X IO^Mq. If the C-image amplification is a bona fide caus- 
tic crossing event, then virial mass measurements underestimate the 
masses of quasars by a factor of several. 



5 DISCUSSION 

5.1 What do caustic crossings really probe for? 

The fine structure of high-amplification events in the light curves 
of microlensed quasars is an effect naturally expected in the frame- 
work of thin accretion disc model. Replacing accretion disc bright- 
ness distribution by some symmetrical structure with brightness 
uniformly decreasing with distance would produce a light curve 
with a single nearly symmetrical peak with the duration time scale 
rd/ve,ff defined by the size of the emitting region (see for e xam- 
ple the analytical solutions given in lSchneider & Weissl ( ll987h . Ap- 
pendix B). Its shape near the maximum is roughly parabolic. In 
this case, the amplification curve maximum is always convex un- 
less the brightness distribution has a non-monotonic or asymmetric 
shape. Interpreting individual humps in the amplification curves as 
individual caustic crossing events leads to unphysically large trans- 
verse velocities several times larger than the upper limits estimated 
in section |3T| 

As it was shown in section 13.21 the only parameter affect- 
ing the amplification curve shape in the simplified disc case is 



X = Td/rin. Disc size scale itself is degenerate with the effective 
transverse velocity. Amplification curves are primarily sensitive to 
the combinations Vin/veff^fJ and rdlv^ffsfl (see section W% 
having the physical meaning of traversal times for the central hole 
in the disc and the disc radial scale itself, respectively. This also 
holds for the general relativistic case as long as the disc is nearly 
face-on. 

Relativistic effects break these degeneracies at inclinations 
i > 60°. Amplification curve properties become sensitive to the 
two angles i and tp, but the dependence on a is weaker and sim- 
ilar light curves emerge for similar values of X if inclination and 
positional angle are the same. 

Interpretation of the three events relies on their consideration 
as caustic-crossi ng events. For QSO J22 3 7+0305, this is j u stified 
by the works of lOil-Merino et all ( l2006h : iKoptelova et^ J2007h 
who performed detailed simulations to check whether artificial am- 
plification maps produce caustic crossings similar (in some sense) 
to the observed peaks in the light curves of individual images. 

For the case of SBS J1520+530, optical depths to microlens- 
ing are relatively small. Central convergence of Kc = 0.5 sets an 
absolute upper limit for it. The optical depth to microlensing is 
smaller than tic unless the stars in the optical path are strongly 
clumped, that is unlikely for an elliptical galaxy (probability of 
accidentally hitting a globular cluster is ~ 10~^). If dark matter 
makes significant contribution to the gravitational field of the lens- 
ing galaxy, optical depths should be decreased proportionally to- 
ward the probable value of ~ 0.1 for the stronger lensed B image. 
Image A is about three times farther from the centre of the lens 
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Figure 11. R-band B/A flux ratio for SBS J1520+530 fitted with a double-fold model (the best-fit curve parameters are given in the text, section l4~Tt . Simplified 
disc single-fold model curve is shown for comparison by a dashed curve. 



that makes it much worse candidate for microlensing. The caus- 
tic net structure in this case is determ ined by small shears distort- 
ing single-lens amplification patterns I IKofman et alJI 19971) . There- 
fore, there is a strong probability for passing a close caustic pair. 
However, in this case the caustics are traversed in opposite direc- 
tions that should produce a pattern completely different from the 
observed amplification curve. 



In this work, we used Kerr-metric geodesies and relativis- 
tic Doppler effects, but the disc model was not fully r elativi stic. 
We used the temperature law of IShakura & SunvaevI ( 119731) in- 
stead of the more a ccurat e relativistic thin disc la w provided by 
iNovikov & Thomel ( Il973h : IPage & Thom3 l ll974l) . Difference in 
the two laws is in the correction factor that decays more smoothly 
inwards and is generally smaller for the fully relativistic thin disc 
case. We tested the difference between the amplification curves cal- 
culated using these two temperature laws and found that for rele- 
vant black hole parameters (M < 10^°Mq and a > -0.9), the 
difference is smaller than one percent of the peak amplification. It 
is more important to consider the effects of non-trivial inner bound- 
ary condition, shock formation and rapid disc tilt change in the in- 
ner parts of the disc and then pass to finer effects like the general 
relativity corrections to the temperature law. 



5.2 Caustic strengths 

For SBS J1520-I-530, amplification curve amplitude is moderate 
and may be characterised by a caustic strength o f k — fii/no ~ 
0.6. Using the estimates made bv lWitt et al.Hl993h . one may expect 
a mean caustic strength of: 



(fc) ~ 0.56 ((m)/r)'''\ 



(20) 



where r ~ 0.1. .0.5 is the optical depth to microlensing, and 
(m) ~ 0.1. .0.5 is the mean microlens mass in Solar units. Note 
that the dependence is very weak both on the optical depth and 
on the mean mass thus making the caustic strength a good con- 
sistency check for the model. In the case of SBS J1520-I-530, ob- 
servational data argue for {m)/T ~ 1, consistent with a moder- 
ately sub-solar mean stellar mass. Low end of the mass function 
of red and brown dwarfs is poorly known and probably variable 
JBastian et aLluOlOi) . There are indications for a low mean stel- 
lar/substellar mas s ((m) ^ 0.1) in the old ste llar population of 
elliptical galaxies jvan Dokkum & Conrovll20I0l) . therefore the op- 
tical depth to microlensing is most likely relatively small, r ~ 0.1. 
If the amplification curve for SBS J1520-I-530 is inteipreted as a 
superposition of two caustic crossings, the two folds have similar 
strengths of 0.7 and 0.6, consistent as well with a stellar population 
rich in low-mass objects. 

For QSO J2237-I-0305 events, the amplitudes are higher. Qual- 
itatively, it was expected, because the lens is a spiral galaxy and the 
number of massive stars among the microlenses should be higher. 
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Figure 12. Amplification curve for the A image of QSO J2237+0305 fitted with a Kerr disc fold traversal model with a = 0.2, Af = 7 X 10^MQ,i = 70°, 
V" = 96°. R band is to the left, V band to the right. Simplified disc best fit is shown for comparison by a dashed curve. Asterisks and circles correspond to 
ISIS and PSF-fitting fluxes, respectively. 



with the mean mass of ~ IMq. Mean caustic strength (k) ~ 2 
(as for the image C event) only if {m)/T ~ IOOMq and hence 
(m) ~ 40Mq (numerical simulations suggest that t he optical 
depth to microlensing is r ~ 0.4 for Einstein's cross tBate et al.l 
UOllh ) that is considerably higher than what one may expect even 
from a young stellar population. In most simulations with r ^ 0.5 
(see fo r example Schechter & Wambsganss (2002): Schechter et al.l 
1 20040 ). predicted amplification factor distributions are broad and 
allow caustics with strengths exceeding the mean value of fc by a 
factor of ~ 2. The reason why our estimates do not contradict the 
theory is selection effect: we consider two brightest and most ev- 
ident high amplification events while the total number of fainter 
events present in the light curves of the four images in the about 
three-year interval of OGLE-II monitoring program is about 10 (if 
the mean distance between fold traversals is about a year). 



5.3 Black hole masses and disc sizes 

The relative size of the inner disc hole and the observed strength 
of relativistic effects favour black hole masses significantly higher 
than virial estimates. For SBS J1520-I-530 and QSO J2237-I-0305, 
the virial masses estimated using the widt h of the CIVA1549 reso- 
nance line are equal to 8.8 x 10^ jPeng et al. 2006) and 9 x IO^Mq 
JYee & De Robertislll99in . correspondingly, with the proposed un- 
certainty of about a factor of 2. Single caustic crossing event fit- 
ting argues for considerably higher masses, (2.. 3) x lO^M© and 
(7..9) X IO^Mq. 

Altering black hole masses by a factor of ~ 3 was proposed 
by [Morgan et alj ( 120 id) to explain the inconsistency between the 



accretion disc sizes estimated using microlensing statistics and the 
accretion disc size predictions following from the standard disc 
theory. Note that in our stud y the degeneracies are different from 
those in [Morgan et alj ( l2010i) . While the Monte Carlo light curve 
analysis is sensitive to the disc size itself (r^ oc {l/rj) ' M^'^), 
the shapes of caustic-traversal events primarily reflect the ratio 
of the disc radial size scale to its inner boundary size (X oc 
{l/riY M~^' ^). To reproduce our results without altering black 
hole masses we should decrease the Eddington ratio by a factor 
of several rather tha n accretion efficiency as it was proposed by 
[Morgan etal] (|20 id). 

To distinguish between higher black hole masses and lower 
I/t] ratios, one may use either photometrical data or different- 
technique accretion disc size estimates. Central black hole mass 
may be estimated using amplification-corrected magnitudes. Ob- 
served flux from a thin disc is calculated by integrating the ob- 
served intensity (indices "obs" and "em" refer here to the radiation 
intensity in the observer frame and in the frame comoving with the 
QSO): 



= Dixll+^y-i COS i X / h^e^ {u/{l+z)) RdR, 

where D — D[z) is angular size distance. This distance scale 
is known for its non-monotonic dependence on redshift. Maximal 
distance of 1.7Gpc is reached at the intermediate redshift of 2 ~ 
1.6 for standard ACDM. Most of the well-studied lensed quasars 
have redshifts close to this value. 
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Figure 13. Best fit for the QSO J2237+0305 image C high amplification event, V band. Simplified disc fit is shown by a dashed line, solid line is the best-fit 
KeiT metric disc fit (a = 0.97, M = 10^°Mq, i = 85° and i/) = 53°). 



For /^, we substitute the Planck law with the temperature law 
determined by the standard accretion disc theory as: 



anc T] 



Here we neglect the correction term that has only small influ- 
ence on the integral flux since the area of the disc affected by the 
term is about (ri„/2.44rd)^ < 0.1 times smaller. For the integra- 
tion limits set to and +oo, the integral is reduced to the following 
JAbramowitz. M. & Stegun. I. A.lll972h : 



/ 



e^ - 1 



= r(8/3)C(8/3) 



Observed flux is finally expressed as: 
© 2009 RAS, MNRAS 000.[Tlfl9l 
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This flux may be used to independently estimate the black hole 
mass. In particular, for the HST F814W filter (we used the calibra- 
tion given in lHoltzman et al.l ( 119950 : this allows to use the magni- 
tudes given by [Morgan et alj ( 1201 Of) in table 1) the mass may be 
estimated photometrically as follows: 



M~2.7x 10^ fj^ 



3/2 



(l + ^)'x 

■^bixio-°-3(^-^9)M0 



(21) 



Following the original work where a similar formula was used 
to estimate the photometrical radius, we denote the magnitude by 
I and use the amplification-corrected values of I(QSO J2237 + 
0305) = 17.9 ± 0.44"" and I{SBS J1520 + 530) = 18.92 ± 
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Figure 14. Mass and l/rj for QSO J2237+0305. Shaded region is tlie band 
restricted by the light curve fitting made in the present worl; (the darkest 
region corresponds to the intersection of the both mass estimates). Vertical 
dashed lines show the range set by the virial mass estimate with the relative 
uncertainty of O.Sdex. Solid lines show the constraints set by equation J211 . 
Disc inclination is set to 60° . 



0.13™. Photometrical masses are closer to the virial mass esti- 
mates. These estimates are degenerate with the l/rj ratio in a way 
different from the results of our light-curve analysis that allows 
to reach consistency if the Eddington ratio is low, l/rj ~ 0.1. 
In figure [141 we show the estimates for l/rj and black hole mass 
made with different methods. Monte-Carlo amplification curve fit- 
ting is sensitive to r^ oc M ' x {l/riY , while the observed 
amplification-corrected flux scales as Fi, oc M*' "^ x {l/rj)^^'^ that 
results in parallel bands in the graph. For QSO J2237+0305, the 
microlensing and photometrical radii are consistent within the un- 
certainties. Our data may be made consistent with photometry and 
virial masses for narrow ranges of parameter values. In particular, 
for QSO J2237+0305, there is a zone where all the three bands in- 
tersect, M ~ (1.6. .1.8) X lO^Mo and l/rj ~ 0.5. .0.6. In the case 
of SBS J1520+530, photometry, virial and light-curve fitting inter- 
sect for M ~ (0.8. .1.1) X lO^Mo and l/rj ^ 0.7..1.3. However , 
in this case the microlensing radius found by [Morgan et al.l ( 120 101) 
is considerably higher. 

Consistency between different methods of radius and mass 
estimates allows to suggest that the microlensing radii found by 
[Morgan et al. (2010) are a subject to some bias that increases in 
some cases the effective radius by a factor of several. This may 
be contamination from a source of much larger angular size (such 
as unresolved starlight or the broad-line region) or some optically- 
thin scatterer. These effects were considered in the original work 
and seem to qualitatively explain the observed discrepancy. 

Since in general there are inconsistencies between the accre- 
tion disc radii obtained by different methods, we refrain from any 
final conclusions upon the masses of the central black holes in QSO 
discs. Still it seems plausible that the apparently large disc sizes 
may be connected to larger black hole masses, smaller Edding- 
ton ratios and high inclinations (the issue of these will be consid- 



ered in detail in the next subsection). We also do not exclude that 
some non-standard accretion regimes may influence the observa- 
tional properties of quasar accretion discs. In particular, low accre- 
tion rates may lead to an optically-thin advection dominated region 
in the inner parts of the disc. The innermost stable orbit radius will 
be over-estimated in this case resulting in black hole masses and 
accretion disc sizes biased toward higher values. We will consider 
the effects of non-standard accretion regimes on microlensing vari- 
ability in subsequent papers. 

Larger black hole masses mean also larger disc sizes. For 
black hole masses M ~ IO'^'^Mq , disc size in the considered wave- 
length range becomes ~ lO^^cm that is comparable to the Einstein 
radius rEin in the case of SBS J1520+530. Straight caustic approx- 
imation is violated in this case because OEin also sets the mean cur- 
vature radius of fold caustics as we ll as the mean d istance between 
individual folds for intermediate r ( [Oaudi & Petters[[2002[) . For the 
case of QSO 12237+0305, Einstein-Chwolson radius is larger be- 
cause of the smaller distance Dl ■ Most of the results of this work 
refer to the inner parts of the disc that makes the eifects of caus- 
tic curvature even smaller. For future more acc urate studies, we 
propo se using the more precise parabolic model dFluke & Webs ten 
[ 19991) . 

5.4 Inclinations of QSO discs 

High inclinations (i > 70°) are required if we interpret the three 
considered high amplification events as caustic crossings. All the 
three events show fine structure near their maxima. It is possible 
to explain this by highly inclined, relativistic discs around black 
holes with masses considerably higher than the virial estimates. 
One amplification curve (QSO J2237+0305, image C) allows to re- 
strict both the Kerr parameter and the mass to a high accuracy. Still 
it should be noted that these estimates are sensitive to the inner 
structure of the disc and therefore model-dependent. High inclina- 
tions are needed because the observational data require asymmetric 
brightness distributions. 

Still, there are strong reasons why QSO discs should not be 
inclined by > 70° : 

• Outer disc rim has a non-negligible thickness and is able 
to obscure the central UV-radiating parts if inclination is i > 
7r/2 — h/R, where h/R ~ 0.02 is the relative disc thickness 
( [Shakura & Sunvaev|[l973l) : this is important if the disc has very 
high inclination i > 85°. 

• A stricter limit is set by dust tori detected spectrally for most 
quasars. Hot dus t emission is observed from QSO J2237+0305 
dAgol et alj [200(1) consistent with a large-scale (R > O.OSpc) 
structure of dust intercepting considerable part of the quasar lu- 
minosity and re-radiating the absorbed luminosity in the infrared 
range. 

• A disc observed edge-on should appear under-luminous. For 
a thin disc, the observed luminosity is cos i times lower than for 
the face-on case, for a concave disc the decrease is stronger. Thus 



the real luminosity should be a factor of cos 



10 higher. In 



the case of QSO 12237+0305, a strong observational evidence for 
small incli nation (i > 60° is ruled out at a 95% confidence level) 
is given in [Poindexter & Kochanela <2010h . 

The strength of the last argument is diminished by the 
high inferred masses of the central black holes. An order of 
magnitude higher black hole mass leads to an order of magni- 
tude higher Eddington luminosity. This makes the observed flux 
consistent with the microlensing disc size estimates. Besides, 
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Figure 15. Flux dependence on inclination for an accretion disc around an 
A/ = IO^Mq, a = 0.6 black hole. Dotted line shows the cosi law for 
compaiison. 
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Figure 16. Simulated amplification curves of an inclined (solid line) and 
tilted (dotted) discs. For both curves, a = 0.9, M = IO^Mq, maximal 



IPoindexter & KochaneM J20ld ) did not consider relativistic effects 
that have the potential to distort the dependence of the observed 
flux on inclination significantly. For the moderate black hole mass 
of IO^Mq and a = 0.6, luminosity dependence on inclination de- 
viates by more than 20% at i > 40° and by more than 50% at 
i > 70° (see figurefTSt. 

An interesting cl ue is that SBS J1520+530 is a BAL (broad ab- 
sorption line) quasar (IChavushvan et al.lll997h . Though BA L QSO 



gh BA L I 

sboooh . ? 



are generally considered as seen nearly edge-on (lElvisbOOOh . some 
of these objects have strongly variable bright radio emission sug- 
esting a c ollimated relativistic outflow observed at low inclination 
Zhou et a l. 2006). It is still an open question how one can coUimate 
a jet in a direction different from the normal to the disc plane. 

A possible solution is that the angular momenta of the accret- 
ing black holes in SBS J1520-I-530 and QSO J2237-I-0305 are mis- 
aligned with the angular momentum of the infalling matter. Since 
there are no reasons for a supermassive black hole to be aligned 
with the accreting matter at large distances, this is probably the 
case. The inner parts of the flow should be tilted and warped. 
Most investigations predict alignment of t he inner parts of the 
flow known as Bardeen-Petterson effect, see lBardeen & PettersonI 
([1975), normally inside the i? ~ 10 ri„ region where the observed 
UV radiation is emitted. Similar misaligrmient may be proposed a 
possible reason for the apparent inconsistency between the inclina- 
tions inferred from radio and optical properties of radio-loud BAL 
quasars. 

In figure [T6] we show the predicted amplification curve for a 
thin accretion disc with zero inclination at infinite distance but tilted 
toward the plane of the black hole inclined by 9bh to the picture 
plane (and to the plane of the disc at large distances). For inclina- 
tion, we use the law 6{r) = di, h ■ exp (—r/ri„ ). Rotation law for 
an inclined disc was taken from lShakural (Il987h . while the temper- 
ature law was assumed standard. Every annulus was considered flat 
that is of course an over-simplification but still serves the main goal 
to show that fold caustic traversal curves are primarily sensitive to 
the innermost parts of the accretion disc. 



inclination is 80°, 
details. 
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- 10000km s-1. See section EU for 



5.5 Influence of disc tilts upon the unification scheme 

In case of random orientation of the black hole spin and the orbital 
angular momentum of the matter entering the disc, an average disc 
will have a tilt of ^ ~ 90° at large distances. Half of the discs 
should be counter-rotating. In case of random misalignment, half 
of the disc planes should have inclinations larger than 60° with 
respect to the black hole rotation frame. A tilted disc around a 
Kerr black hole is usually proposed to gradually change its inclina- 
tion with decreasing radius approachin g alignment at the innermos t 
stable orbit Linear analysis made by IZhuravlev & Ivanovl ( 1201 ih 
for a thin disc around a nearly-Schwarzschild black hole confirms 
this conclusion only for the case of considerably large viscosity 
(see the original work for details). The authors also show stabil- 
ity of counter-rotating discs around slowly rotat ing black holes. On 
the other hand, numerical simul ations made by [Fragile & AnninosI 
( l2005h : lDexter & Fragile ( 1201 ll) suggest that the disc evolution to- 
ward alignment should take place in its innermost parts, at several 
Schwarzschild radii. 

Tilted discs may become a serious obstacle for formation and 
observability of relativistic jets. Once formed, relativistic outflows 
will be stopped by the higher-pressure material of the disc or the 
dust torus. If jet formation and propagation into the rarefied galac- 
tic medium becomes impossible in case of strongly inclined discs, 
one may try to explain the observed radio brightness dichotomy 
of QSO by different regi mes of alignment, as it was proposed by 
[Lawrence & Elvisl ( 120101) for Seyferts and radiogalaxies. Presence 
of a thick torus with 8t ~ 30° leads to about 75% probability that 
the jet will be stopped either by the torus or by the tilted disc itself. 
The real fra ction of radio-weak sources is higher, about 90% (see 
for example lHewett et alj ( 1200 ll) ). 

Importance of disc tilts for jet formation may also have con- 
nection to the formation of the "localised" X-ray radi ation found in 
some QSO. In particular, in the two recent works bv lZimmeret al.l 
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( I2OIII) and IChen et all J201lh . the properties of the X-ray emit- 
ting region in QSO J2237-I-0305 are recovered using microlens- 
ing curves. In IChen et alj ( 1201 ih . it is concluded that this region 
is smaller and more variable than any existing extended jet or 
corona model can explain. Among other possibilities, the "failed 
jet" model by Ghise llini et al.l (2004) is put forward. If the jet pos- 
sesses enough energy to escape the gravitational well but is stopped 
by the disc matter, the source of the X-ray emission may be con- 
nected to the hot spot in the disc where its energy is dissipated. 
Relativistic termination shock remains invisible but the bow shock 
where the ram pressure is balanced by the pressure of the disc 
should make important contribution to the observed X-ray emis- 
sion. And the spatial size of the X-ray bright spot may be small if 
the jet is sufficiently collimated. Stopped jet model may be checked 
by considering correlations between different properties of X-ray 
and radio emission from quasars. 
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6 CONCLUSIONS 

We show that general relativity effects have considerable influence 
on the amplification curves of microlensed quasars. This effect is 
more profound if the disc is strongly inclined. Having better data 
on hand (smaller observational errors, larger homogeneity and bet- 
ter temporal coverage) may allow to resolve the structure of high- 
amplification events with a better precision and to use them for 
accurately measuring disc tilts and black hole rotation parameters 
of lensed quasars. 

Some implications may be made even now. In particular, the 
discs in both objects, SBS J1520-I-530 and QSO J2237-I-0305, are 
seen at high inclinations. Apparent contradiction with the high ob- 
served fluxes and lack of strong absorption signatures in the spec- 
tra may be resolved if one considers a disc having low inclination 
(i < 60°) at tens of GM/(? but is seen nearly edge-on in its inner 
parts. Similar picture is expected if the initial angular momentum 
of the accreted matter is strongly misaligned with that of the black 
hole. 

High probability to effectively measure a high inclination may 
be qualitatively understood as a very strong tilt and warp near the 
last stable orbit. In this case, there should be, at a high probabil- 
ity, some radius at which the disc lies perpendicular to the picture 
frame. Contribution of this annulus should be enhanced by Doppler 
boost. In other words, if a disc has a range of inclinations due to 
tilts and waips, higher inclinations will make larger contribution to 
microlensing amplification curves. 

In the analysed data we do not find any challenges for the stan- 
dard accretion disc model. But they bring our attention toward the 
role of disc tilts that was never considered in the framework of QSO 
microlensing. Besides, some of accretion disc properties are better 
understood if we alter the accretion disc parameters such as Ed- 
dington ratio and accretion efficiency. The apparently high black 
hole masses and disc sizes required by microlensing effects may 
be connected to the low accretion efficiency and contamination by 
larger angular scales. 



NOTE ADDED IN PROOF: 



We became aware of the papers by iBogdanov & CherepashchukI 
( I2OO2I 120040 where amplification curves considered in this work 
were analysed. Authors restore one-dimensional brightness distri- 
butions and propose to connect their shapes to relativistic effects. 
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